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A REVIEW OF RECENT RESULTS ON SPREAD F THEORY 


Introduction 


discovered some forty years ago by Booker and Wells (1938). Until 


recently, there has been only much statistical data concerning 


Spread F. However, in the last few years with advances in radar 


and numerical simulation techniques a clearer picture of the funda 


mental plasma instability mechanisms causing equatorial Spread F 


phenomena has been evolving. It should be emphasized that this paper 


will deal only with equatorial Spread F theory and, in addition, dwell 


only on those theories using plasma mechanisms as a basis. For 


experimental results, the works of Farley et al. (1970), Dyson et al 


(1974), Kelley et al. (1976), Woodman and La Hoz (1976), Morse et al 


(1977), and McClure et al. (1977) should serve as good references for 


In speaking about Spread F one should at least show its basic 


Spread F emerged from the results of ionosonde traces such as those 


exhibited in Fig. 1. The multiple traces are caused by magnetic field 


aligned irregularities. If there were no irregularities, a single 


crisp trace rising to the right would be exhibited on the ionogram 


At this point we want to remember the basic equatorial geometry 


Figure 2 exhibits the basic equatorial nighttime ionospheric F region 


geometry, i.e., the geometry under which equatorial Spread F occurs 
Note: Manuscript submitted November 16, 1978. 










N (y) represents the background electron density as a function of 


altitude (y). Gravity, £, points down, the ambient magnetic field, B, 
is horizontal (pointing north) and k represents a horizontal perturba¬ 
tion (in the westward direction). The maximum in the electron density 
profile is the F peak. The underside of the profile steepens at night 
due to chemical recombination effects and electrodynamic forces. The 
E region has been severely reduced by chemical recombination and plays 
a negligible role. To a plasma physicist this geometry is a classi¬ 
cal flute mode geometry and one might expect this equatorial iono¬ 
spheric geometry to be unstable to a variety of plasma instabilities. 
However, until recently, one of the basic difficulties was getting the 
unstable irregularities to the topside (i.e., above the F peak) when 
they were initiated on the bottomside. One must remember that the 
experimental evidence has exhibited both top and bottomside 
irregularities. 

At this point a general brief review will be given of equatorial 
Spread F (ESF) theories. First, we will discuss the linear theories. 
Dungey (1956) was the first to suggest that ESF was initiated on the 
bottomside by a Rayleigh-Taylor instability. Dagg (1957) suggested 
that the ESF phenomena was due to E to F region coupling, i.e., 
irregularities in the E region coupled up to the F region. In 1959, 
Martyn (1959) was the first to suggest that ESF was a manifestation of 
the E X B gradient drift instability. Calvert (1963) proposed that 
the downward motion of the neutral atmosphere at night was responsible 
for ESF. This mechanism is essentially equivalent to the E X B 








instability because of the relative motion between ions and neutrals 
in determining the instability. All of these previously invoked 
linear instability mechanisms could only explain the formation of 
bottomside irregularities. The collisional Rayleigh-Taylor instability 
with field line averaging was also proposed (B alsley et al. , 1972; 
Haerendel , 1974) as a linear instability mechanism. By averaging 
(integrating) the density along the magnetic field, the total electron 
content profile becomes steeper on the bottomside and its peak is 
raised in altitude with respect to the local electron density peak. 

This would allow the linear mechanism to operate to slightly higher 
altitudes ('v 100 km greater) , but still would not explain the 
existence of irregularities above this "new peak". Hudson and 
Kennel (1975) pointed out the importance of the collisional drift 
mode in ESF in the wavelength regime 30m - 100m. This mode could be 
excited on both the top and bottomside but still would not explain the 
longer wavelengths. In their paper, finite larmor radius (FLR) 
corrections were also applied to the collisionless and collisional 
Rayleigh-Taylor instability. 

Several nonlinear theories have been invoked to explain the 
different ESF observations. For example, Hudson et al. (1973) sug¬ 
gested that the very smallest scale (^.10m) irregularities (e.g., 
those seen by radar coherent backscatter) were due to a two step 
process. In this prescription a longer wavelength instability sets 
up the driving conditions for the shorter wavelengths to become 
unstable. This is similar in spirit to the successful two step 






theory ( Sudan et al. , 1973) proposed for Type II equatorial E region 
electrojet irregularities. Haerendel (1974) suggested that the range 
of wavelengths (many kilometers down to centimeters) exhibited by ESF 
phenomena was due to a multi-step process. This scenario is as 
follows: (i) the collisional Rayleigh-Taylor (R-T) instability with 

horizontal wavevectors is driven by gravity and the background, zero 
order electron density gradient scale length on the bottomside; then 
(ii) the E X B gradient drift instability with vertical wavevectors 
arises due to the horizontal density, large amplitude variations set 
up by the collisional R-T instability; then (iii) the inertia 
(collisionless) dominated R-T instability arises; and finally (iv) 

kinetic drift waves grow upon these irregularities after they reach 

-2 

large amplitude. Chaturvedi and Kaw (1976) interpreted the k 
measured power spectrum of the ESF plasma density irregularities in 
terms of a two step theory. In this theory longer wavelength R-T 

modes couple to kinetic collisional drift waves in such a manner that 

-2 

the mode coupling results in the observed k spectrum. 

I 

A major breakthrough was made by Scannapieco and Ossakow (1976) 
who performed a nonlinear numerical simulation of the collisional R-T 
instability for ESF geometry. The simulation results showed that the 
collisional R-T instability generated irregularities and bubbles 
(plasma density depletions) on the bottomside of the F region which 
subsequently rose beyond the F peak by nonlinear polarization induced 
E^ X B forces. This was the first theoretical result to explain how 
long wavelength irregularities could appear on both the bottomside and 
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topside of the F region. The bubble phenomena was in accord with the 
recent observations ( Kelley et al. t 1976; McClure et al. , 1977, 

Woodman and La Hoz , 1976) of plasma density depletions. An analytical 
nonlinear mode-mode coupling theory for the coherent development of 
the collisional R-T instability was performed by Chaturvedi and 

Ossakow (1977). This theory suggested that vertical modes would be 

-2 

dominant and result in a k power spectrum. Hudson (1978) extended 
the previous results to the collisionless R-T regime and reached 
similar conclusions. Analytical models for the rise of collisional 
and collisionless R-T ESF bubbles, in analogy with fluid bubbles, 
was presented by Ott (1978). At the same time, Ossakow and Chaturvedi 
(1978) presented analytical models for the rise of collisional R-T ESF 

I 

bubbles within the context of the electrical analogy with barium 
clouds. 

Costa and Kelley (1978a,b) suggested that coherent steepened 

-2 

structures and not turbulences would give a k power spectrum. More¬ 
over, these sharp gradients could cause small scale sizes (^ 20m) by 
collisionless low frequency (much less than the ion gyrofrequency, ft^) 
kinetic drift waves via a two step process. Their analysis was a 
linear one carried out on a nonlinear state, i.e., one achieves the 
steepened gradients by nonlinear processes and then one performs 
linear theory on this state. Kelley and Ott (1978) suggested that the 
ESF bubbles, in the collisionless R-T regime, generate a wake with 
vortices. They then applied two dimensional fluid turbulence theory 
to the model. This resulted in the development of turbulence at 
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shorter and longer wavelengths than the bubble size. This in turn led 
to a prediction of k ^ for the power spectrum (which does not appear 
to be in agreement with existing experimental observations) in the 
range I. g 1 < k < L Q \ where L g is the stirring (bubble) size and L^ 
is a dissipation length cutoff. In a continuation of the numerical 
simulation work, Ossakow et al. (1978) showed a more rapid ESF develop¬ 
ment and higher bubble rise velocities resulting from sharper bottom- 
side background electron density gradients and higher altitudes of the 
F peak. In Huba et al. (1978) very small scale (wavelengths lm and 
36cm) irregularities are reported. A two step process, utilizing high 
frequency ,(^(K) kinetic drift cyclotron or lower hybrid drift 
instabilities, is invoked to explain them. Linear theory for these 
instabilities, in the ESF environment, was performed on the nonlinear 
state with encouraging results. 

The above introductory remarks glaringly point out that much work 
and significant progress in the theoretical area of equatorial Spread 
F has been accomplished in the past few years (indeed just look at 
the number of publications during 1978 alone). Notwithstanding the 
recent successes, much work still needs to be done. Indeed, the 
theoretical and numerical simulation efforts in ESF are continuing 
along a hot and heavy path. Section II of this paper presents out¬ 
lines of some of the theoretical efforts briefly mentioned in the 
preceding paragraphs. Given the length limitations, it would be 
exceedingly difficult to outline all of the theoretical works 
mentioned or even to present all of the details of a few works. 
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Section II, hopefully, will wetten the reader's appetite to read the 
referenced works. Section III presents a summary concerning ESF 
theory. 



J = ne (V. - V.) 
— —i —e 


(5) 


In the above equations the subscript a denotes species (e is electron, 
i is ion), n is density, V is velocity, P is the production, v R is the 
chemical recombination rate, T is temperature, V is the gradient 








operator, i;. the ambient magnetic field (taken to be uniform) , 
e is the electronic charge, c is the speed of light, m is mass, g. is 
gravity, v. is the ion-neutral collision frequency, J is current, and 
the electrostatic approximation has been made where E = —V4>. 

Equation (1) is the continuity equation, (2) and (3) are the electron 
and ion momentum equations, respectively, (4) is the divergence of 
the current and (5) is the current equation. What we have in mind is 
to apply the set of equations (1) - (5) to the two dimensions perpen¬ 
dicular to at the geomagnetic equator, making various approxima¬ 
tions. 

Assuming a harmonic perturbation dependence of the form 

exp-i (k • x - oot) , where x denotes perpendicular to B (k is hori- 
_ — x —o — 

zontal) and linearizing equations (1) - (5) we obtain for the linear 
growth rate 

r 2 Vn o 15 1 

Y= -v. + (v. - 4 g • -) I/2 - v_ ... 

L in m * n J' R (6) 


oj r a) + iy 
r 

which reduces to 



L| 


2 

v . >>4 g/i 

in 

v*\ << 4 g/L 
in 

I Vn I 


O 


n 


o 


(7a) 

(7b) 


B 







in the collision dominated and inertia dominated regimes, for 
kL>> 1. Equations (7a) and (7b) represent the R-T instability, 
including recombination damping, in the respective regimes. (It 
should be noted that this only represents instability on the bottom- 
side of the F region where the first term in (7a) and (7b) is 
positive.) 

b. 2D Computer Simulation Results . In this section we will 
outline the two dimensional ( J- B^) computer simulation results 
( Scannapieco and Ossakow , 1976; Ossakow et al. , 1978). The simula¬ 
tions follow the nonlinear evolution of the collisional R-T 
instability; consequently, in eqn. (3) inertial terms (i.e., the left 
hand side of the eqn.) are neglected. Furthermore, one takes for the 
F region /fi^« 1 (Q^ = eB Q /m^c) and the quasi-neutrality assumption 

is made, i.e., n a; n. a n. Equations (1) - (5) then become 
e l 



V • J = 0 , J = ne (V. - V ) (11) 

— — —l -e 

where one has taken T g = T^ = 0 for simplicity (see Ossakow et al. , 
1978), B^ = B q z and the subscript o refers to equilibrium quantities. 
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and breaking the potential into an equilibrium and a perturbed 


quantity, 


linear growth rate 


This clearly shows that only the bottomside of the F region where 


Vn is positive can be linearly unstable (and only if the first term 


Equations (14) and (15) were solved numerically using a 


vertical mesh spacing of Ay= 2 km and total y extent of 200km, and 


an east-west, horizontal mesh spacing Ax - 200m and a total horizontal 


and v as a function of 






altitude were utilized. The system was initialized with a 
perturbation of a few percent in the horizontal (x), east-west 
direction with a wavelength , v3km and the evolution in time of (14) 
and (15) was followed for different background electron density 
profiles. Figure 3 shows the results for a background electron 
density, n Q , profile with an F peak at 354km and a minimum bottomside 
background electron density gradient scale length, L '''10km. In this 
case at t = 4000 sec a bubble (plasma density depletion) is clearly 
forming and beginning to rise in the central portion of the mesh 
(note n = n + n^). The isodensity contours are such that the maximum 
absolute value of the enhancement or depletion is in the center and 
the contours decrease (in absolute percentage) as one goes toward the 
outer contours. At t = 4000 sec the maximum depletion within the 
bubble is 54% and the maximum enhancement over the mesh is 84%. At 
t = 8000 sec we notice that the bubble has reached the altitude of 
the F peak, with the innermost contour of the rising bubble repre- 

4 

senting a 41% depletion. At t = 10 sec the main bubble is clearly 
through the F peak with an innermost depletion contour of 41%. How¬ 
ever, in the ionosphere below the bubble near x = 0 there is a 71% 
depletion contour, similarly in the wings near |x| = 4km. The inner¬ 
most enhancement contour, at this time, represents a 236% enhancement 
with a maximum inside this contour of 294%. Note that the top of the 
main bubble is at an altitude of 375km while the bottom trail of the 

4 

bubble is at an altitude of 270km. Between t = 8000 a: d 10 sec the 
bubble has risen '''24km which represents a rise velocity '''12m/sec. 











Also note that the bubble is ''-1km wide. Figure 4 depicts contours 

4 

of constant induced potential ^ at t - 10 sec. This shows that 
the more isolated part of the high altitude bubble depicted in Fig. 3 
is acted on by an induced electric field which points from west to 
east and is dipolar in nature. This causes the bubble to rise with a 
(-c/B q ) V x z velocity. However, the lower portion of the mesh is 
acted on by an induced electric field •■'hich points from east to west. 
This field is much weaker than the induced field acting on the 
isolated portion of the central bubble. The lower altitude electric 
field causes the enhancements and depletions to move downward. Thus, 
the lower altitude portion of the central bubble becomes captured by 
the enhancement?. 

Figure 5 presents for comparison a case in which the background 
electron density profile (shape) was kept the same, but the entire 
profile moved up in altitude so that the F peak was at 434km. All 
other parameters are the same as in Figure 3, except and v R are 
taken for the altitude range 332km to 532km (those used in this 
simulation). One can immediately note the more rapid time evolution 
of the Spread F process with respect to that presented in Fig. 3. 

At t = 700 sec, a rising bubble with an innermost contour of 79% and 
a maximum depletion inside this contour of 84% were noted. At 
t = 1000 sec the bubble has reached the peak and at this time the 
innermost depletion contour is 85%. The long trail associated with 
the high altitude bubble has a 100km extension to lower altitudes. 

At t = 1400 sec the top of the main bubble is at an altitude ^SOOtan 






and has a long trail connecting to an altitude of 357km. There is a 
maximum 70% depletion within the innermost contour of the high 
altitude bubble. Between t = 10 3 and 1400 sec the top part of the 
bubble rose ' x '65km and this represents a rise velocity 'V'lBOm/sec. 
Potential contour results for this simulation show similar patterns to 
those exhibited in Fig. 4. Naturally, the induced electric fields 
causing the bubble to rise in the present case are stronger. This 
spread in bubble rise velocities has been observed by AE satellite 
data ( McClure et al. , 1977). 

Other numerical simulations in this series have been performed 
and the paper by Ossakow et al. (1978) should be consulted for more 
details. The basic conclusions reached from these simulations are as 
follows: (i) the collisional R-T instability causes linear growth on 

the bottomside of the equatorial Spread F region; (ii) plasma density 
depletions (bubbles) steepen on their top and nonlinearly rise to the 
topside ionosphere, beyond the F peak, by polarization (induced) 

E X B forces; and (iii) high altitude of the F peak, small bottomside 
background electron density gradient scale lengths, and large percent¬ 
age depletions yield large vertical bubble rise velocities, with the 
first two conditions favoring collisional R-T linear growth (instabil¬ 
ity) . In addition, large spatial bubbles with similar rise velocities 
to those presented here, but with almost 100% depletions, have been 
produced by numerical simulations ( Zalesak et al. , 1978). In these 
cases the horizontal mesh covered 200km in extent (intermesh spacing 
Ax = 5km) and a long wavelength (''-75km) initial perturbation was 
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(18) 


~ = A. , tsin (k v-u>t) cosk x + A_ sin 2k x 
n 1,1 y x 2,0 x 

o 

is chosen. In this analysis x is vertical (altitude) and y horizontal 
(east-west). This is the way it appears in the reference. For the 
convenience of the reader, in referring to the original work, we have 
kept the coordinate system of the reference. 


The coupled mode equations for the amplitudes A^ ^ and A^ g 


become 

3A 


1,1 


3 t 


= Y 


1 1 ^ 1 1 ” ^ ^ 1 1 A 2 0 


(19) 


3A 


3t 


2,0 , A a „ 2 

Y 2,0 A 2,0 2 A l,l 


( 20 ) 


where the coupling coefficient a = k k * g/k^ v. , and u> = ui + iy. 

x y in r ' 

It should be noted that the linear growth rate y^ 0 ne 9 at i- v e» i.e., 
the purely vertical mode Q is linearly damped. Also y^ ^ is posi¬ 
tive and represents linear growth of the mode A 1 ^ In the saturated 
steady state one has SA^/St = 3A 2 Q /3t = 0 and from (19) and (20) 
this results in 


' 1,1 


2,0 


1,1 


2 k L 
x 

>5 


( 21 ) 


;/k 2 Vin\ /- 2Y 2,0 a \ 

W^7) 


( 22 ) 






(i.e., the power spectrum would k_ 


coherent nonlinear evolution where the linearly damped mode A. 


generated nonlinearly by the linearly unstable mode A. 


monic generation (A contains k and A_ n contains 2k ). More 


detailed information can be found in Chaturvedi and Ossakow (1977) 


d. Analytical Models for ESF Bubbles. First we will discuss 


collisional R-T bubbles (Ossakow and Chaturvedi, 1978) in which the 


electrical analogy with plasma density enhancement (e.g., barium 


clouds) was utilized and results obtained for general 2D < J § 0 ) 


bubble shapes. Here eqns. (14) and (15) are utilized with the 


further simplifying assumptions of neglecting recombination chemistry 


0) and the explicit altitude dependence of the ion-neutral 


The following set of equations result 


where an ambient horizontal electric field E has been included to 


show more generality. Equation (24) can be thought of, in general 


electric field, E*. The plasma density depletion (bubble) is 


analogous to the case of a cavity immersed in a dielectric with a 











uniform electric field (E*). In general (23) and (24) do not admit of 


a two dimensional f jB ), steady state solution. However, for the case 


For a constant density depletion with an elliptical shape one has 


where H is the Heaviside function (note the geometry is taken such 


that x is vertical and y horizontal). For this elliptical shape 


the solution to (24) , neglecting E , is 


Using -(c/B ) x z this further results in the nonlinear vertical 


bubble rise velocity, V , given by 


Limiting cases of (29) for sheet, cylindrical and slab bubble 











geometries are given by 


9/v. 

in 


- 6n/n Q ) 1 , 

b « a 

(30a) 

- 6n/n Q ) 1 , 

b = a 

(30b) 
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b << a 

(30c) 


For typical values of v in as a function of altitude. Fig. 6 exhibits 

rise velocities given by (30b) for various percentage depletions. 

Table 1 shows the rise velocities, in units of g/v. , for various 

in 

shapes. (Note that the linear case comes from linearizing (23) and 

(24)). All of the geometry results can be expressed in a concise 
formula 

9 ..Si 


v 

B v. n 
in o 


(31) 


where f(5n/n Q ) is an increasing function of the percentage 
depletion 6n/n Q . Basically, the results predict that high altitudes 
and/or large percentage depletions yield high vertical rise velocities 
for the bubbles (in agreement with experimental observations). 

Collisional and collisionless (inertia dominated) two dimensional 
cylindrical R-T bubbles have been studied analytically by Ott (1978). 
This study is based on the analogy with fluid dynamic flows and brings 
forth some of the work done on two-dimensional fluids. This study 
begins with the basic equations (1) - (3), considers two dimensions 
(j3q) , sets n e % n^ss n, and makes the assumption that 


a i ” St - 


7, v 


in 


(31) 


To lowest order, using (31), one obtains from (3) a lqwest order ion 
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velocity (with z = B /1B |), 


(Note that inks units are being used here to coincide with the units 


used by Ott (1978)). Quasi-neutrality, i.e., V • J = 0, with the 


assumption of two dimensionality implies that J can be specified in 


terms of a single scalar potential function, i|>, such that 


Using eqn. (2) for the electron velocity, and a next order ion 


velocity equation obtained by putting (32) into (3), the following 


ion velocity equation is obtained 


(o) „ , v> (o) 

= - Vp - nm.v. V. 

l in—l 


with p = p + p. + B <J». Using (32) and the assumption if either 


isothermal or adiabatic ions one has 


Using in the ion continuity equation (1) with v R = P = 0, 


eqn. (1) becomes 











Equations (34) - (36) form a complete set of equations 

(o) % 

sufficient to determine the unknown quantities , n, and p. In 

the limit of v. -*• o, eqns. (34) - (36) are identical to those of 
an ideal incompressible fluid. At this point the philosophy taken 
is that there is much to be learned concerning bubbles in the ESF 
ionosphere from the extensive studies of bubbles in fluids. One 
then uses a stream function £ such that 


V. = z x V £ 


(37) 


After a series of manipulations one can obtain an equation for the 
stream function. In the collisional and collisionless case one 
finds that for certain values of thq bubble rise velocity £ will 
satisfy the equation, for a cylindrical shape (see Ott, 1978 for 



The results of this study by Ott (1978) predicts the bubbles 


to be cylindrical (circular cap at top) in two dimensions. The 


bubble rise velocity in the collision-dominated regime, for a 100% 


depletion, is given by 


which is altitude dependent. This result is the same as that pre 


1, i.e., a 100% depletion. In the 


dieted by (30b) for 6n/n 


inertia (collisionless) dominated regime, for a 100% depletion, the 


bubble rise velocity is given by 







where R is the radius of curvature at the top of the bubble so that 
(39) is size dependent. 

e. Analytical Work on the Very Small Scale (^10m) 

Irregularities . The work in this section represents essentially 
multilinear calculations using kinetic theory resulting in plasma 
kinetic drift wave modes. It is multilinear because it depends on 
say a two-step process whereby linear theory is performed on the non¬ 
linear state. The driving density gradients, in these calculations, 
are thought to arise through a primary instability, driven by the zero 
order background ionospheric equatorial F region electron density 
gradient, achieving a large amplitude state. The zero order iono¬ 
spheric electron density gradient is of larger scale size than the 
primary instability electron density gradient scale size, which would 
be of the order of the instability wavelength. Because the calcula¬ 
tions are kinetic, they employ particle distribution functions. 

Kinetic drift waves have been investigated for laboratory plasma 
fusion conditions for over twenty years, so a well developed formalism 
exists. 

Before proceeding with specific calculations let us present some 
general concepts regarding kinetic drift waves which will be useful 
for both types of calculations presented in this section. The basic 
geometry is such that 

B = B z (40a) 

— —o 

n = n (x) (40b) 

o o 
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(similar to the zero order background equatorial ionospheric 


l 





geometry) where here n Q (x) arises due to the primary instability. The 
orthogonal coordinate system is completed by the ion and electron 
diamagnetic velocities being along the y axis, 

^d - (v di - V de> y 

V J . = (V. 2 /2Q.) (31n n /3x) 

di 11 o 

V, = - (V 2 /2Q ) (31n n /3x) 

de e e o 

i 

where V. = (2T. /m. ) and the larmor radius is defined by 

i e ife i e 

r * V /ft . The linear analysis is then performed with this in mind, 
a a a 1 

First the low frequency (u) « ft^) collisionless drift wave 
calculations of Costa and Kelley (1978a,b), for kr^ 'v 1 (r^ is the ion 
gyroradius), will be presented. These calculations were meant to pro¬ 
vide a basis within which to try and account for the 3 meter radar 
backscatter observations of Jicamarca (see Woodman and La Hoz , 1976). 

A linear kinetic dispersion relation is derived using perturbations of 
the form exp i |^k • x - wt j, where k = k^ z + k y and w = + iy. 

Furthermore, assumptions are made such that w « ft., V. <<|w/ki, |« V 

1 i || e 

and u> > v. , v. , v.., which are the ion-neutral, ion-electron, and 
in le n ' 

ion-ion collision frequencies, respectively. This analysis depends on 
having a k ^ , a component of the wavevector along the ambient geo¬ 

magnetic field. Figure 7 shows some growth rate results from these 
calculations. The growth rate is in units of the ion thermal velocity 


(41b) 

(41c) 
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divided by the electron density gradient scale length 

(L = n^ (dn Q /dx) . Figure 8 depicts some measured inverse electron 

density gradient scale lengths during ESF and maximum growth rates as 

a function of these inverse scale lengths. The basic results of these 

calculations show maximum growth for k r. a 1.5, i.e., X s: 20m (for 

-L 1 ± 

typical ESF parameters) with growth rates ^ 1 sec 1 . For more 
detailed analysis concerning these low frequency drift waves applied 
to ESF see Costa and Kelley (1978a,b). 

Now we present the high frequency drift wave analysis of Huba 
et al. (1978). In this reference radar backscatter observed 
irregularities with wavelengths of 1 meter and 36 cm at Kwajalein dur¬ 
ing ESF conditions are shown. In an effort to explain these very short 
wavelength irregularities high frequency (to ^ ) drift waves were 

analyzed for ESF conditions. These waves are the so-called drift 
cyclotron (DC) or lower hybrid drift (LHD) instabilities with maximum 



growth rates for k r ^ 1 (and k. 


these instabilities. The parameter determining which instability 


Utilizing the linear dispersion relation for high frequency drift 


L/r. < {1/21) (m./m ) 


where i is the harmonic number (u> 


plasma this requires the electron density gradient scale length 










L < 340m, which is satisfied (see Fig. 8). Growth rates for 
these instabilities are given by 


Y 25 (m /iik) £ 0^ 


( 44 ) 


2 7 

However, the condition C, << 1 implies that (k r.) n << 2 x 10 and 
£ ± 1 

for k r 'v 1 this means n « 10 3 cm 3 , which is quite restrictive, 
j. e 

Longer wavelengths, i.e., smaller values of k r ^ would raise the 
density restriction somewhat, but still be restrictive for ESF 
conditions. 

For C f ^ 1, the lower hybrid drift instability is operative and 
there is no threshhold condition on L. Basically, the collisions 
which increase C f destroy the ion gyroresonances needed for the DC 
instability to operate. The real and imaginary part of the fre¬ 
quency for the LHD instability are given by 


*5 ,2 1 

u) r./m.\ /r. \ / m. \ 

_J1 ^ _L(_L ) , ( -± ) 

«i L fi. \L ) \m e J 


(45) 



:n the collisionless limit it should be noted that the DC instability 


transforms into the LHD instability for high enough ion diamagnetic 


velocities such that 


which for ESF conditions implies that L ^ 30m. Figures 9 and 10 show 


some typical results from the analysis of the high frequency drift 








The results of this analysis predicts that the lower hybrid 


drift instability is dominant for most typical ESF ionospheric para¬ 
meters. Also, maximum growth of the instability occurs for 

k r <\i 1 (X <v 21 cm) , although good growth rates can occur for 
J. e x 

X 'x, 1 m. Finally, from this instability, large growth rates 


x 

(y fi.) resulting in growth times, t = y less than a second can 
occur. For more details of this work see Huba et al. (1978). 


III. Summary 

Although much progress has been made in the theoretical efforts 
directed toward the equatorial Spread F ionosphere, especially in the 
past three years, more has to be done. Also the burden cannot be 
placed on the theoretician alone. Correlative measurements have to 
be made prior to and during ESF conditions. It is not sufficient to 
make a single measurement with one instrument and then expect a 
complete theoretical description of ESF. One needs to know the state 
of the ionosphere with respect to driving parameters such as back¬ 
ground electron density profiles (bottomside electron density 
gradient scale lengths and height of the F peak), d.c. electric 
fields, neutral winds, and ionic mass composition in order to build a 
predictive model. In addition, in order to compare results from the 
predictive model, measurements of the in situ fluctuating component 
of the electric field and plasma density have to be made, as well as 
radar backscatter measurements of the very small scale irregularities 
(^10m) and ground measurements of satellite signal propagation 
amplitude and phase scintillations. 










Some of the remaining theoretical problems could be listed as 
follows. (1) What are the effects of changing initial conditions, 
including ion inertia, and including neutral winds in the numerical 
simulations? (2) An analytical description of many bubbles rising. 

(3) How do bubbles decay and what role does diffusion, etc. play? 

(4) How does a turbulent development occur? (5) What are the effects 
of other regions of the ionosphere (e.g., E region) on ESF? (6) A 
more complete study of collisional effects on drift waves is needed. 
What determines when the small scale irregularities should occur (a 
more quantitative description)? (7) A determination of the nonlinear 
saturation of small scale irregularities (instabilities) is needed. 

(8) What are the effects of ? Indeed points (6) and (7) are tied 
to the more general question of what is the relation between the very 

small scale (^10m) irregularities (which radar backscatter observes) 
and the longer wavelength fluid type (e.g., R-T) irregularities (which 
are primarily responsible for Spread F seen on ionograms)? I am sure 
that some more questions and points could be raised. However, the 
above list should keep the theoreticians busy for a reasonable time. 
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Fig. 3 - Contour plots of constant nj/n Q for the simulation with an F 
peak at 354km at t * 2000, 4000, 8000, and 10,000 sec. The small dashed 
contours with a plus sign inside and the solid contours with a minus sign 
inside indicate enhancement and depletions over the ambient electron 
number density. The large dashed curve depicts the ambient electron 
number density (values on upper horizontal axis), n Q , as a function of 
altitude. The vertical y axis represents altitude, the lower horizontal 
x axis is east-west range, and the ambient magnetic field is along the 
z axis, out of the figure. Taken from Ossakow et al. (1978) 


















Fig. 4 - Contours of constant induced (polarization) potential, 0]. 
over the mesh (see Fig. 3) at t = 10,000 sec. Plus and minus denote 
positive and negative values, with values decreasing in magnitude as 
one goes from the innermost to the outermost contours. The large 
dashed curve is n Q . Taken from Ossakow et al. (1978) 
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Fig. 7 - Growth rates as a function of b for the low frequency collisionless drift waves 
(Pi is the ion gyroradius). The curve Y << oo r is obtained from a small growth rate approx' 
imation expression. Taken from Costa and Kelley (1978b) 
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INVERSE GRADIENT SCALE LENGTH (dlnn/dx) OF PLASMA DENSITY, meters 

Fig. 8 - Observed Inverse gradient scale lengths (1/n) (dn/dx) as a function of altitude 
during the Natal rocket flight. Also indicated by the straight line is the maximum growth 
ratecf low frequency drift waves as a function of (1/n) (dn/dx). Note that p. is the ion 
gyroradius. Taken from Costa and Kelley (1978a) 1 
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Fig. 10 - Maximum growth rate as a function of ion diamagnetic 
drift velocity for the DC instability (solid line) and the LHD 
instability (dashed line). The other parameters are the same 
as in Fig. 9. Note that the DC instability goes into the IHD 
instbillty for Vjj/Vi ?> 0.11. Taken from Huba et al. (1978) 
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(■>PAC£ D£ c ENSc SYSTEMS 
OlCY ATTN SZJ 


safoiA laboratories 

P. 0. BOX 5800 
ALBUQUERQUE, 8’l15 

OlCY ATTN DOC CON FOR J. MARTIN 

OlCY ATTN DOC CON FOR W. BROWN 

OlCY ATTN DOC CON FOR A. ThORNBROUG n 

OlCY ATTN DOC CON FOR T. WRIGHT 

OlCY ATTN DOC CON FOR D. DAhLGREN 

OlCY ATTN Da CON FOR 3141 

OlCY ATTN Da CON FOR SPACE PROJECT DIV 







SANOIA LABORATORIES 
LIVERMORE LABORATORY 
P. 0. dOX %69 
LIVERMORE, CA 94550 

OiCY ATTN OOC CON FOR d. MURPHEY 
01CY ATTN OOC CON FOR T. COCK 

OFFICE OF MILITARY APPLICATION 
department of energy 
WASHINGTON, D.C. 20545 

OiCY ATTN OOC CON FOR 0. GALE 


OTHER GOVERN NT 


CENTRAL INTELLIGENCE AGENCY 
ATTN RO/SI, RM 5G48, HQ BLOG 
WASHINGTON, D.C. 20505 

OICY ATTN OSI/PSIO RM 5F 19 

DEPARTMENT OF COMMERCE 
NATIONAL bureau of standards 
WASHINGTON, D.C. 20254 

CALL CORRES: ATTN SEC OFFICER FOR) 
OICY ATTN R. .MOORE 

OEPARTI'CNT OF TRANSPORTATION 
OFFICE OF THE SECRETARY 
TAD-44.1, ROOM 10402-3 
400 7Th STREET, S.W. 

WASHINGTON, D.C. 20590 
OICY ATTN R. LEWIS 
OICY ATTN R. DOHERTY 

INSTITUTE FOR TELECOM SCIENCES 
NATIONAL TELECOMMUNICATIONS & 1^0 AOMIN 
BOULOER, CO 80305 

OiCY ATTN A. JEAN (UNCLASS ONLY) 
OICY ATTN W. UTLAUT 
OICY ATTN D. CROMBIE 
OICY ATTN L. BERRY 

NATIONAL OCEANIC & ATMOSPHERIC AOMIN 
cNVIRPMCNTAL RESEARCH LABORATORIES 
DEPARTMENT OF COWERCE 
BOULDER, CO 80302 

OICY ATTN R. GRUBB 

OiCY ATTN AERONOMY LAB G. REID 


DEPARTMENT OF DEFENSE CONTRACTORS 


AEROSPACE CORPORATION 
P. 0. BOX 92957 
LOS ANGELES, CA 90009 

OICY ATTN I. GARFUNCEL 
OICY ATTN T. SALMI 
OICY ATTN V. JOSEPHSON 
OICY ATTN S. BOWER 
OICY ATTN N. STOCKWELL 
OICY ATTN 0. OLSEN 
OICY ATTN J. CARTER 
OICY ATTN F. MORSE 
OICY ATTN SMFA FOR PWW 

ANALYTICAL SYSTEM6 ENGINEERING CORP 
5 OLD CONCORD ROAD 
BURLINGTON, W 01803 

OICY ATTN RADIO SCIENCES 

BERKELEY RESEARCH ASSOCIATES, INC. 
P. 0. BOX 983 
BERKELEY, CA 94701 

OICY ATTN J. WORKMAN 


BOEING COMPANY, THE 
P. 0. BOX 3707 
SEATTLE, wA 98124 

OICY ATTN G. KEISTER 
OiCY ATTN D. MURRAY 
OICY ATTN G. HALL 
OICY ATTN J. KENNEY 

CALIFORNIA AT SAN DIEGO, UNlV OF 
IPAPS, B-019 
LA JOLLA, CA 92093 

OICY ATTN HENRY G. BOOCER 

BROWN ENGINEERING COH 5 ANY, INC. 

CUMMINGS RESEARCH PARK 
HUNTSVILLE, Al 35807 

OICY ATTN ROMEO A. DELIBERIS 

CHARLES STARK DRAPER LABORATORY, INC. 

555 TECHNOLOGY SQUARE 
CAMBRIOGE, MA 02159 

OICY ATTN 0. B. COX 
OICY ATTN J. P. GILMORE 

COrf»UT£R SCIENCES CORPORATION 
6565 ARLINGTON 3LVD 
FALLS CHURCn, VA 22046 
OICY ATTN H. BLAMC 
OICY ATTN JOHN SPOOR 
OICY ATTN C. NAIL 

COMSAT LABORATORIES 
linthicum road 
CLARKSBURG, MD 20734 
OICY ATTN G. HYDE 

CORNELL UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 

ITHACA, NY 14850 

OICY ATTN 0. T. FARLEY JR 

ELECTROSPACE SYSTEMS, INC. 

BOX 1359 

RICHARDSON, TX 75060 

OICY ATTN H. LOG STUN 

OiCY ATTN SECURITY (PAUL PHILLIPS) 

ESL INC. 

495 JAVA DRIVE 
SUNNYVALE, CA 94066 

OiCY ATTN J. ROBERTS 
OICY ATTN JAMES MARSHALL 
OICY ATTN C. W. PRETTIE 

FORD AEROSPACE & CCmJNlCAT IONS CORP 
3939 FABIAN WAY 
PALO ALTO, CA 94303 

OICY ATTN J. T. MATTINGLEY 

GENERAL ELECTRIC COMPANY 
SPACE DIVISION 
valley FORGE SPACE CENTER 
GOOOARD BLVO KING OF PRUSSIA 
P. 0. BOX 8555 
PHILADELPHIA, PA 19101 

OICY ATTN M. H. BORTNER SPACE SCI -Afl 

GENERAL ELECTRIC CO*»ANY 
P. 0. BOX 1122 
SYRACUSE, NY 13201 

OICY ATTN F. RE1BERT 

general electric co*»any 

TE^O-CENTER FOR AOVANCED STUDIES 
816 STATE STREET (P.O. DRAWER QQ) 

SANTA BARBARA, CA 93102 
OiCY ATTN DASIAC 
OICY ATTN DON CHANDLER 
OiCY ATTN TOM BARRETT 
OICY ATTN TIM STEPHANS 
OICY ATTN WARREN S. <WP 
OICY ATTN WILLIAM MCNAMARA 
OICY ATTN 3. GAMBILL 
31CT ATTN MACK STANTON 
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ieKRAL ELECTRIC TECH SERVICES CO. 
hMES 

COURT STREET 
Syracuse, ny 13201 

01CY ATTN G. HILLMAN 


<AMAN SCIENCES CORP 
P. 0. BOX 7463 
COLORADO SPRINGS, CO 80933 
OiCY ATTN T. MEAGHER 


L1Mt A6IT COR? 

10**53 ROSELLE 
SAN 01 EGO, CA 92121 

OICY ATTN IRWIN JACOBS 


GEKRAL RESEARCH CORPORATION 
SANTA BARBARA DIVISION 
P. 0. BOX 6770 
SANTA BARBARA, CA 93111 

OICY ATTN JOHN ISE JR 
OICY ATTN JOEL GAR0ARI NO 


LOWELL RSCH FOUNDATION, UNIVERSITY OF 
450 AIKEN STREET 
LOWELL, MA 01854 

OICY ATTN <. BIBL 


GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AT. 99701 

CALL CLASS ATTN: SECURITY OFFICER) 
OICY ATTN T. N. DAVlS CuNCL ONLY) 

OICY ATTN NEAL BROWN (UNCL ONLY) 

OICY ATTN TECHNICAL LIBRARY 


M.l.T. LINCOLN LABORATORY 
P. 0. BOX 73 
LEXINGTON, HA 02173 

OICY ATTN DAVID M. TOWLE 
OICY ATTN P. WALDRON 
OICY ATTN L. LOUGHLlN 
OICY ATTN 0. CLARK 


GTE SYLVANIA, INC. 

ELECTRONICS SYSTEMS GRP-EASTERN DIV 
77 A STREET 
NEEDHAM, MA 02194 

OICY ATTN MARSHAL CROSS 


MARTIN .MARIETTA CORP 
ORLANOO DIVISION 
P. 0. BOX 5837 
ORLAhOO, FL 32805 

OICY ATTN R. HEFFNER 


ILLINOIS, UNIVERSITY OF 

departknt of electrical engikering 

JRBANA, IL 61803 

OiCY ATTN K. YEH 


MCDOWELL DOUGLAS CORPORATION 
5301 BOLSA AVENUE 
HUNTINGTON BEACH, CA 92647 
OiCY ATTN N. HARRIS 
OICY ATTN J. MOULE 
OICY ATTN GEORGE MR02 
OICY ATTN W. OLSON 
OICY ATTN R. *. HALPRIN 
OICY ATTN TECHNICAL LIBRARY SERVICES 


ILLINOIS, UNIVERSITY OF 
107 COBLE HALL 
801 S. WRIGHT STREET 
URBANA, IL 60680 

(ALL CORRES ATTN SECURITY SUPERVISOR FOR) 
OICY ATTN K. YEH 


INSTITUTE FOR DEFENSE ANALYSES 
400 ARMY-NAVY DRIVE 
ARLINGTON, VA 22202 

OICY ATTN J. M. AEIN 
OICY ATTN ERNEST BAUER 
OiCY ATTN HANS WOLFHARD 
OICY ATTN JOEL BENGSTON 


MISSION RESEARCH CORPORATION 
735 STATE STREET 
SANTA BARBARA, CA 93101 
OICY ATTN P. FISCHER 
OICY ATTN W. F. CRcVlER 
OICY ATTN STEVEN L. GUTSCME 
OICY ATTN D. SAPPEDIELD 
OICY ATTN R. BOGUSCH 
OICY ATTN R. HENDRICK 
OICY ATTN RALPH <ILB 
OICY ATTN DAVE SOWLE 
OICY ATTN F. FAJEN 
OICY ATTN M. SCME10E 
OICY ATTN CONRAD L. LONGMIRE 
OICY ATTN WARREN A. SCHLUETER 


MSS, INC. 

2 ALFRED CIRCLE 
BEDFORD, MA 01730 

OICY ATTN DONALD HANSEN 


INTL TEL 4 TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLET, NJ 07110 

OICY ATTN TECHNICAL LIBRARY 


MITRE CORPORATION, THE 
P. O. BOX 208 
BEDFORD, MA 017 30 

OICY ATTN JOW .MORGANSTERN 
OICY ATTN G. HARDING 
OICY ATTN C. E. CALLAHAN 


JAYCOR 

1401 CAMINO DEL MAR 
DEL MAR, CA 92014 

OICY ATTN S. R. GOLDMAN 


MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADI SON BLVD 
MCLEAN, VA 22101 

OICY ATTN W. MALL 
OICY ATTN W. FOSTER 


JOWS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOWS HOPKINS ROAD 
LAUREL, *> 20810 

OICY ATTN DOCJMENT LIBRARIAN 
OICY ATTN THOMAS POTEMRA 
OICY ATTN JOW DASSOULAS 


PACIFIC-SIERRA RESEARCH CORP 
1456 CLOVERFI ELD BLVD. 

SANTA MONICA, CA 90404 

OICY ATTN E. C. FIELD JR 


LOOKED MISSILES 4 SPACE CO INC 
P. 0. BOX 504 
SUWYVALE, CA 94088 

OICY ATTN DEPT 60-12 
OICY ATTN 0. R. CHURCHILL 


PEMVSYLVANIA STATE UNIVERSITY 
IONOSPHERE RESEARCH LAB 
318 ELECTRICAL ENGINEERING EAST 
UNIVERSITY PARK, PA 16802 

(NO CLASSIFIED TO THIS ADORESS) 

OICY ATTN IONOSPHERIC RESEARCH LAB 


LDOCKED MISSILES AW SPACE CO INC 
3251 HANOVER STREET 
PALO ALTO, CA 94304 

OICY ATTN MARTIN WALT DEPT 52-10 
OICY ATTN RICHARD G. JOWSON DEPT 52-12 
OICY ATTN W. L. I*KDF DEPT 52-12 








PHOTOMETRICS, INC. 

442 MARRETT ROAD 
LEXINGTON, *A 02173 

OICY ATTN IRVING L. KOFSKY 

PHYSICAL DYNAMICS INC. 

P. 0. SOX 3027 
BELLEVUE, WA 93009 

01CY ATTN E. J. FREMOUW 

PHYSICAL DYNAMICS INC. 

P. 0. BOX 1069 
BERKELEY, CA 94701 

01CY ATTN A. THOMPSON 

R S 0 ASSOCIATES 

p. 0. BOX 9693 

MARINA DEL REY, CA 90291 

01CY ATTN FORREST GILMORE 
01CY ATTN BRYAN GA88ARD 
01CY ATTN WILLIAM B. WRIGHT JR 
01CY ATTN ROBERT F. LELEVlER 
01CY ATTN WILLIAM J. KAR2AS 
01CY ATTN M. ORY 
01CY ATTN C . MACDONALD 
01CY ATTN R. TuRCO 

RAfO CORPORATION, Th£ 

1700 MAIN STREET 
SANTA MONICA, CA 90406 

01CY ATTN CULLEN CRAIN 
01CY ATTN ED BEDROZIAN 

RIVERSIDE RESEARCH INSTITUTE 
80 WEST END AVENUE 
NEW YORK, NY 10023 

01CY ATTN VINCE TRAPANI 

SCIENCE applications, INC. 

P. 0. BOX 2351 
LA JOLLA, CA 92038 

01CY ATTN LEWIS M. LlNSCN 
OICY ATTN OANIEL A. HAM.IN 
01CY ATTN D. SACHS 
OICY ATTN E. A. STRAKER 
OICY ATTN CURTIS A. SMITH 
OICY ATTN JACK MCDOUGALL 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 

OICY ATTN BARBARA ADAMS 

SCIENCE APPLICATIONS, INC. 

HUNTSVILLE DIVISION 
2109 W. CLINTON AVENUE 
SUITE 700 

HUNTSVILLE, AL 35805 

OICY ATTN DALE H. DIvIS 


SCIENCE APPLICATIONS, INCORPORATED 
8400 WESTPARK DRIVE 
MCLEAN, VA 22101 

OICY ATTN J. COOCAYf* 

SCIENCE APPLICATIONS, INC. 

80 MISSION DRIVE 
Pleasanton, ca 94566 
OICY ATTN SZ 


SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MEN.0 PARK, CA 94025 

OICY ATTN DONALD NEILSON 
OICY ATTN ALAN BURNS 
OICY ATTN G. SMITH 
OICY ATTN L. L. COBB 
OICY ATTN DAVID A. JOHNSON 
OICY ATTN WALTER G. CHESNUT 
OICY ATTN CHARLES L. RINO 
OICY ATTN WALTER JAYE 
OICY ATTN M. BARON 
OICY ATTN RAY L. LEAOA8RAND 
OICY ATTN G. CARPENTER 
OICY ATTN G. PRICE 
OICY ATTN J. PETERSON 
OICY ATTN R. HAKE, JR. 

OICY ATTN V. GONZALES 
OICY ATTN D. MCDANIEL 

TECHNOLOGY international CORP 
75 WIGGINS AVENJE 
BEDFORD, MA 01730 

OICY ATTN W. P. BOQUIST 

TRW DEFENSE t SPACE SYS GROUP 

ONE SPACE PARK 

REDONOO BEACH, CA 90278 

OICY ATTN R. K. PLEBUCH 
OICY ATTN S. ALTSCHULER 
OICY ATTN D. DEE 

V1SIDYNE, INC. 

19 THIRD AVENUE 

north west industrial park 

BURLINGTON, MA 01803 

OICY ATTN CHARLES HUMPHREY 
OICY ATTN J. W. CARPENTER 
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IONOSPHERIC MODELING DISTRIBUTION ^1ST 
UNCLASSIFIED ONLY 


PLEASE DISTRIBUTE ONE COPY TO EACH OF Th£ FOLLOWING PEOPLE: 


ADVANCED RESEARCH PROJECTS AGENCY (ARPA) 
STRATEGIC TECHNOLOGY OFFICE 
ARLINGTON, VIRGINIA 


HARVARD UNIVERSITY 
HARVARD SQUARE 
CAMBRIDGE, MASS. 02138 


CAPT. DONALD M. LEV IKt 


M. B. MCELROY 
R. LINDZEN 


NAVAL RESEARCH LABORATORY 

WASHINGTON, D.C. 20375 PEl'WSYLVANI A STATE UNIVERSITY 

UNIVERSITY PARK, PENNSYLVANIA 16802 


DR. P. MANGE 
DR. R. MEIER 

DR. E. SZUSZCZEWICZ - COOE 7127 

DR. TIMOTHY COFFEY - COOE 6700 (20 COPIES) 

DR. S. OSSAKOw - COOE 6780 (100 COPIES) 

DR. d. GOODMAN - COOE 7560 


SCIENCE APPLICATIONS, INC. 
1250 PROSPECT PtAZA 
LA JOLLA, CALIFORNIA 92037 

DR. O. A. HAM.IN 
DR. L. LINSON 
DR. D. SACHS 


DR. J. S. NISBET 
DR. P. R. ROHRBAUGh 
DR. D. E. 3ARAN 
DR. L. A. CARPENTER 
DR. M. LEE 
DR. R. DIVANY 
DR. P. BENNETT 
DR. E. KLEVANS 


UNIVERSITY OF CALIFORNIA, lOS ANGElES 

405 HIlLGARD AvENUc 

LOS ANGELES, CALIFORNIA 90024 

DR. F. V. CORON IT| 

DR. C. KEMVEL 


DIRECTOR OF SPACE AND ENVIRCNhCNTAu LABORATORY 
NOAA 

BOULDER, COLORADO 80302 


UNIVERSITY OF CALIFORNIA, BERKELEY 
BERKELEY, CALIFORNIA 94720 


dr. *. Glen* jean 

DR. G. W. ADAMS 
DR. D. N. ANOERSON 
DR. <. DAVIES 
DR. R. F. DOMVELLY 


A. F. GEOPHYSICS LABORATORY 
L. G. HANSOM FIELD 
BEDFORD, MASS. 01730 

DR. T. ELKINS 
DR. w. Sw!^R 
''OS. R. SAGALYN 
DR. J. M. FORBES 
DR. T. J. KENEShEA 
DR. J. AARONS 


OFFICE OF NAVAL RESEARCH 
800 north quincy street 
ARLINGTON, VIRGINIA 22217 

DR. H. MULLANEY 


DR. M. HUDSON 


UTAH STATE UNIVERSITY 
4Th N. ANO 8Th streets 
LOGAN, UTAH 84322 

DR. P. M. BANKS 
DR. R. HARRIS 
DP. PETERSON 
DR. R. MEGILL 
DR. K. BAKER 

CORNELL UNIVERSITY 
IThACA, NEW YORK 14850 

dr. w. e. Swartz 

DR. R. SUDAN 
DR. D. FARLEY 
DR. M. KELLEY 
DR. E. OTT 


NASA 

GOOCARD SPACE FlIC^I CENTER 
GREENBElT, MARYLAM) 20771 


COMMANDER 

NAVAL electronics laboratory center 
SAN DIEGO, CALIFORNIA 92152 

DR. M. BLEIwEISS 
DR. 1. ROThKJLLER 
OR. V. HILDEBRAND 
MR. R. ROSE 


•j. S. army ABERDEEN RESEARCH ANO DEVELOPMENT CENTER 
ballistic research laboratory 
ABERDEEN, MARYLAND 

DR. J. HEI.MfRL 


DR. S. CHANORA 
DR. K. MAEDO 


PRINCETON UNIVERSITY 

plasma physics laboratory 

PRINCETON, NEw JERSEY 08540 

DR. F. PERKINS 
DR. E. FRYMAN 


INSTITUTE FOR DEFENSE ANAcYSIS 
400 ARMY/NAVY DRIVE 
ARLINGTON, VIRGINIA 22202 

DR. E. BAUER 


COMMA »OER 

NAVAL AIR SYSTEMS COWMANO 
DEPARTMENT OF THE NAVY 
WASHINGTON, O.C. 20360 

DR. T. CZuBA 


UNIVERSITY OF PITTSBURGH 
PITTSBURGH, PA. 15213 

DR. N. ZABUSKY 
DR. M. 91JNOI 


49 











